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Massive black holes are believed to reside at the centres of most galax- 
ies. They can become detectable by accretion of matter, either contin- 
uously from a large gas reservoir or impulsively from the tidal disrup- 
tion of a passing star, and conversion of the gravitational energy of the 
infalling matter to light. Continuous accretion drives Active Galactic 
Nuclei (AGN), which are known to be variable but have never been ob- 
served to turn on or off. Tidal disruption of stars by dormant massive 
black holes has been inferred indirectly but the onset of a tidal disrup- 
tion event has never been observed. Here we report the first discovery 
of the onset of a relativistic accretion-powered jet in the new extra- 



galactic transient, Swift J164449.3+573451. The behaviour of this new 
source differs from both theoretical models of tidal disruption events 
and observations of the jet-dominated AGN known as blazars. These 
differences may stem from transient effects associated with the onset 
of a powerful jet. Such an event in the massive black hole at the centre 
of our Milky Way galaxy could strongly ionize the upper atmosphere 
of the Earth, if beamed towards us. 

Swift J164449.3+573451 (hereafter Sw J1644+57) was discovered by 
the SwiyPBurst Alert Telescope^ (BAT, 15-150 keV) in a 20-minute image 
of the sky that began at 12:57:45 UT on 28 March 201 P. The Swift X-Ray 
Telescope^ (XRT) began observing the field at 13:22: 19.8 UT and found a 
bright, uncatalogued X-ray source. The most accurate XRT positiorP^is 
RA (J2000): 16 h 44 m 49.92 s , Dec (J2000): +57° 35' 00.6", with a 90% 
confidence error circle radius of 1.4 arcseconds. A precise localization by 
the EVLA found a variable source at the centre of a galaxy within the XRT 
error circle^. Optical spectroscopy of this soured^ measured a redshift of 
0.354, corresponding to a luminosity distanc43of 5.8 x 10 27 cm. 

Subsequent analysis of BAT data taken before the on-board trigger 
shows that the outburst was first detected on 25 March 2011, with peak 
BAT count rates of ~ 0.02 cts cm -2 s -1 during the days just before the 
on-board trigger (Supplementary Figure 1 Beginning on 28 March the 
BAT data show multiple flares peaking at up to 0.09 cts cm -2 s _1 (about 
20% of the count rate from the Crab nebula). The source flux then dropped 
dramatically, with an average count rate of 0.0020±0.0005 cts cm -2 s _1 
between 2 April 201 1 and 18 April 201 1. 

The flares seen by BAT are tracked in the 0.3-10 keV band by the XRT 
(Figure 1, Supplementary Figures 1 and 2). Following nearly 3 days of in- 
tense flaring with peak fluxesjover 5 x 10~ 9 erg cm~ 2 s _1 , Sw J1644+57 
decayed over several days to a flux of about 8 x 10 -12 erg cm -2 s _1 , then 

(Figure 2). Over the next two 



rose rapidly to ~ 8 x 10 erg cm s 
weeks the X-ray flux decayed gradually to ~ 1.5 x 10 



erg cm 



'assuming Ho = 71 km s L , Q m = 0.27, and Qa = 0.73. 

"During the brightest parts of the outburst, Sw J1644+57 had power law photon indices 
ranging from 1.3 to 1.8 in the BAT bandpass; the corresponding flux conversion factor is ~ 
10 7 erg/count. See Supplementary Table 1. 
3 All XRT fluxes are for the 0.3-10 keV band. 
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with occasional sharp dips to < 3 x 10 -11 erg cm -2 s _1 . No statistically 
significant periodic or quasi-periodic signals were found in the XRT data 
(see Supplementary Methods). 

The X-ray spectral shape varies strongly with X-ray flux in the first 
few days after the BAT trigger (Supplementary Figures 3 and 4). A simple 
power law can fit the spectra, but the fits are improved with one of several 
more complex models (see Supplementary Methods for details). Substan- 
tial absorption in the host galaxy is required, with Nh ~ 2 x 10 22 cm -2 
(Supplementary Figure 5). 

We obtained photometry in the uvw2, uvm2, uvwl, u, b, v, R, J, H, 
and Ks bands, as well as broad-band "white" with the Swift UVOT, TNG, 
BOAO, LOAO, and Maidanak Observatory telescopes (see Supplementary 
Methods for details). The heavily obscured optical counterpart of the tran- 
sient X-ray source is not detected in optical or UV bands, where the host 
galaxy dominates the light, but is detected strongly in the near-infrared. 

Upper limits from the Fermi Large Area TelescoptP (LAT) and from 
VERITAS^prowide constraints in the GeV and TeV band, respectively (see 
Supplementary Methods for details). 

Archival X-ray and 7-ray data show no evidence of Sw J1644+57 be- 
fore 25 March 2011. X-ray flux upper limits from observations between 
1990 and March 24 are orders of magnitude lower than the peak X-ray 
fluxes observed by BAT and XRT, and are significantly lower than the 
lowest flux recorded by XRT since Swift observations began on 28 March 
2011 (see details in Supplementary Methods). The source has brightened 
by at least 4 orders of magnitude compared with some earlier deep obser- 
vations of this field, and we consider that the source was truly dormant for 
decades before this outburst began. 

The abrupt appearance of a bright new persistent extragalactic X-ray 
source is unprecedented. Gamma-ray bursts reach similar fluxes, but fade 
rapidly much more rapidly than Sw J1644+57. Blazars (AGN with ener- 
getic jets pointed toward Earth) are known to be highly variable, but are 
persistent and have never been observed to turn off or to turn on, or display 
this level of variability. Supemovae in distant galaxies are rarely detected 
at X-ray wavelengths. 

The isotropic X-ray luminosit)|3 of Sw J1644+57 (corrected for ab- 
sorption in the source) ranges from 10 45 — 4 x 10 48 erg s _1 , with an 
integrated X-ray power 4 (over the 23 days following the BAT trigger) of 
~ 3 x 10 53 erg. If Sw J1644+57 radiates isotropically, its total energy 
production is significantly larger than the electromagetic energy released 
by a supernova explosion. The persistence of the X-ray emission implies 
that this object is not powered by an impulsive event such as an exploding 
star; instead, it requires a continuous production of prodigious amounts of 
energy. Only gravitational energy released by accretion of matter in the 
vicinity of a black hole can produce such power levels over such a long 
time period. 

The observed minimum 3er X-ray variability time scale of 
Sw J1644+57, St hs ~ 100 s, constrains the size of the black hole un- 
der the assumption that the central engine dominates the var iabilitjB For 
a Schwarzschild black hole with mass Mhh and radius r s , the minimum 
variability time scale in its rest frame is <5£ m i n ~ r s /c ~ 10.0(Me) s, 
where M 6 = (M bh /W S M @ ). At z = 0.354, this gives 



7.4 x 10° 



5to 



100 s 



Mr. 



(1) 



Much smaller masses are unlikely. The black hole mass can also be con- 
strained through the Mbh — ibuigc relation, which gives an upper limit 
of 2 x 10 7 Mq. We conclude that the black hole mass is in the range 
1 < Ma < 20 (see Supplementary Discussion). 

The Eddington luminosit)0 of this black hole is Z/Edd = 1.3 x 
10 44 ergs _1 M6 . The observed peak luminosity is super-Eddington 

4 In the XRT band, which is 0.4-13.5 keV in the Sw J1644+57 rest frame. 
Variability can also be produced in the jet, as is seen in blazars.. 

6 The luminosity at which the gravitational force on infalling gas is balanced by outward 
radiation pressure. 



and requires a strongly anisotropic radiation pattern with a relativistic 
jet pointed towards us. A relativistic jet is also required by the radio 
transienP, which we interpret as an external shock in the gas surround- 
ing Sw J1644+57. 

The spectral energy distribution (SED) produced from our data (Fig- 
ure 3) shows that the broad-band energy spectrum is dominated by the 
X-ray band, which accounts for 20% (50%) of the total bolometric en- 
ergy output in the low (high/flaring) X-ray state 4 . The SED places strong 
constraints on the emission mechanism. The IR-to-X-ray slope, nearly as 
hard as v 1 ' 3 , can only be accommodated within synchrotron theory in the 
absence of low-energy electrons, which would severely overproduce the 
optical flux. This requires a particle-starved source, implying that the jet 
is magnetically-dominated (see Supplementary Discussion for details, and 
discussion of other models). At high energies, the upper limits from Fermi 
and VERITAS require that the self-Compton flux from the X-ray peak is 
absorbed by the 7-7 — > e process. The bulk Lorentz factor in the X-ray 
emitting region must therefore be T < 20; we assume a value of ~ 10. 

This luminous, accretion-powered, relativistic jet can be powered by 
the tidal disruption of a star or by the onset of AGN activity from a pre- 
viously dormant massive black hole. The Swift BAT, with a field of view 
of ~ 47r/7 sr, has detected one such event in ~ 7 years at a flux that 
would have been detectable to z ~ 0.8. The all-sky rate of Sw J1644+57- 
like events is therefore ~ 1 yr -1 , with a 90% confidence interval 
of 0.08 — 3.9 yr - . The rate of tidal disruption events is thought to 
be ~ 10~ 5 yr -1 galaxy -1 on both observational and theoretical 
grounds. Given the galaxy number density n ga i ~ 0.01 Mpc~ 3 , the total 
number of galaxies in the co-moving volume within z < 0.8 can be esti- 
mated as A?gai(z < 0.8) ~ 10 9 . This gives a total tidal disruption event 
rate of R T UAn(z < 0.8) ~ 10 4 yr" 1 . If these act like AGN, 10% of them 
contain relativistic ("radio-loud") jets. The fraction of these jets that point 
towards us is 

or 



max 



2T 



2 • 



where V is the bulk Lorentz factor and 8j is the opening angle of the 
jet. A bulk Lorentz factor of V ~ 10 — 20 or a jet opening angle of 
6j ~ 5° would provide the observed event rate. Our observations are 
consistent with theoretical predictions of the formation of a a low-density, 
magnetically-dominated jet during the super-Eddington phase of a tidal 
disruption event around a 10 6 — 10 7 Mq black hole^. In the tidal disrup- 
tion interpretation, we expect the source to begin a slow decay once the 
material in the accretion disk is exhausted and the luminosity begins to 
track the fall-back timescale (t -5//3 ) of the disrupted stellar material onto 
the black hole^. This decay may have begun a few days after the BAT 
trigger (Figure 2). 

A dormant massive black hole could also become visible through the 
onset of long-lived AGN activity. The typical life time of an AGN is 
tAGN ~ 3 x 10 7 yfQ In steady-state, the rate at which AGN turn on 
is therefore ~ i^GN ( DU|; see Supplementary Discussion). Within the 
z < 0.8 volume, the event rate of the onset of a new radio-loud AGN 
would be i?AGN,tot = 0.1JV ga i(z < 0.8)/£agn ~ 3 yr -1 . A moderate 
Lorentz factor of V ~ 3 or an opening angle of Oj ~ 13° would have 
an event rate of ~ 0.1 yr - , consistent with the observations. The onset 
could be initiated by an accretion disk in a rotating galactic bulge, popu- 
lated by gravitationally captured mass lost from stars in the inner galaxy 
through supernovae or stellar winds. The inner radius of this disk will de- 
crease slowly under the effect of a viscosity, with a final infall timescale 
of td ~ 2 x 10 3 (a/0.1) _1 M6 s, resulting in a rapid onset of electro- 
magnetic radiation. In this case one expects an extremely long-lived event 
(~ £agn). Long-term monitoring of Sw J 1644+57 will be required to 
distinguish between these scenarios for the accretion. 

We can consider Sw J1644+57 to be an analog of Sgr A*, the massive 
black hole in t he ce ntre of our Milky Way galaxy. Sgr A* has a mass of 
~4x 10 6 MeP^I, and yet has a quiescent soft X-ray luminosity^EH of 



only 4 x 10 32 erg s -1 , which would be undetectable at z=0.354. If a jet 
similar to Sw J1644+57 turned on in Sgr A* and was pointed at Earth, the 
peak X-ray flux at the Earth could reach ~ 200 erg cm -2 s _1 . The X- 
ray flux would correspond to a Y200-class solar flare, > 40 times brighter 
than the strongest X-ray flare ever detected from the SurP, and the fluence 
(integrated flux) could exceed 10 4 J m~ 2 . Such an event would produce 
strong ionization effects on the Earth's atmosphere, and the total fluence is 
> 10% of that thought necessary to produce mass extinctions in the case 
of a 10 s duration 7-ray bursP^. Fortunately, the probability of a Sgr A* 
jet being pointed at Earth is low. 
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Figure 1 Swift BAT (top) and XRT (bottom) light curves of the early 
flares from Sw J1 644+57. The horizontal axis is time since the BAT 
trigger on 28 March 201 1 . The BAT and XRT count rates track each 
other, with episodes of flaring in the first few days after the BAT trigger 
on 28 March 201 1 . Data gaps are caused by times when the source 
was not being observed. Error bars are one standard deviation. 
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Figure 2 Swift XRT light curve of Sw J1 644+57 for the first 3 weeks 
of observations. The X-ray events were summed into time bins contain- 
ing 200 counts per bin and count rates were calculated for each time bin. 
Because of the strong spectral evolution of this source and its large in- 
trinsic absorption, light curves in physical units cannot be produced with 
constant count rate-to-flux conversion factors. We use time-dependent 
conversion factors to convert count rates to fluxes. Spectral fits were per- 
formed on time slices of data that typically contained about 3000 counts 
(see Supplementary Methods). Count rate-to-flux conversion factors 
were calculated for each spectrum, using absorption-corrected fluxes 
in the rest frame 0.4-13.5 keV band. These time-dependent flux con- 
version factors were then interpolated to obtain fluxes for each time bin. 
The right-hand axis gives the conversion to luminosity of the source, as- 
suming isotropic radiation. The dashed line shows a t ~ 5/3 curve (where 
t is time since the BAT trigger). Error bars are one standard deviation. 
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Figure 3 | The Spectral Energy Distribution (SED) of Sw J1 644+57. 

The green data points show the transient brightness a few days after 
the BAT trigger; cyan data points show the measured SED at 4.5 days; 
black data points are XRT data between 6.5 and 9.5 days after the BAT 
trigger. Near-infrared (NIR) were obtained with the TNG and BOAO tele- 
scopes in the Ks band 48 hours after the trigger, and in the J, H, and Ks 
band fluxes 4.5 days after trigger. The NIR flux has been dereddened 
with A v = 4.5. Details of the observations and data reduction are 
presented in Supplementary Methods. In the X-ray band we show the 
spectrum at the peak of the bright flares (31 hours after BAT trigger; 
XRT and BAT) and in the low flux state (~ 4.5 days after, XRT only), 
together with the spectrum in the intermediate persistent flux state that 
began a week after the BAT trigger (time integrated between day 6.5 
and 9.5 from trigger, XRT only). The X-ray data have been corrected for 
absorption with a constant N H = 2 x 10 22 cm~ 2 . Upper limits from 
the Fermi LAT at 2 x 10 23 Hz and from VERlTA^at 10 26 Hz are also 
shown. The red curve shows the model discussed in the text, which is a 
blazar jet model^fit to our SED. The dominant emission mechanism is 
synchrotron radiation peaking in the X-ray band. On the low frequency 
side, the hard slope between the NIR and X-ray bands requires sup- 
pression of low-energy electrons, which would otherwise overproduce 
the NIR flux. On the high frequency side, the LAT and VERITAS upper 
limits require that the self-Compton component is suppressed by 7-7 
pair production. The model includes a disk/corona component from the 
accretion disk (black dotted curve), but the flux is dominated at all fre- 
quencies by the synchrotron component from the jet. The blue curve 
shows the corresponding model in the low X-ray flux state. The kink in 
the X-ray spectrum suggests that a possible additional component may 
be required; it would have to be very narrow, and its origin is unclear. 
The model parameters for this fit are given in the Supplementary Dis- 
cussion, and suggest a magnetically-dominated jet, which is unusual for 
AGN. This might be explained by the unusually high accretion rate and 
small black hole mass (compared to typical AGN) implied by our data. 
In the SI we discuss two alternative models: one with a magnetically- 
dominated jet but a more AGN-like parameter set, and one at the oppo- 
site extreme of a particle-dominated jet. Both of these alternative models 
seem contrived, and we favor the model presented here. All error bars 
are one standard deviation. 
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Discovery of the Onset of Rapid Accretion by a Dormant 
Massive Black Hole: Supplementary Information 



Supplementary Methods 

Swift Discovery and broad-band observations 

BAT data The Swift 25 Burst Alert Telescope 26 (BAT, 15-150 keV) triggered on a new uncatalogued source at Tq = 
12:57:45 UT 28 March 201 \Zk The trigger was assumed to be a y-ray burst (GRB) and was named GRB 1 10328A, 
following standard nomenclature. However, this trigger was followed by three additional BAT triggers at increasing 
count rates over the next two days 28 (see Supplementary Table [TJ, making it clear that this object was not, in fact, a 
GRB. (Following the fourth trigger, the on-board BAT source catalog was adjusted to prevent further triggers.) The 
X-ray counterpart was renamed Swift J164449. 3+573451, which is now the preferred name of the source 29 . We will 
refer to Swift J 164449. 3+573451 as Sw J 1644+57 for the remainder of this paper. 

Post-facto examination of pre-trigger data indicate that the outburst was first detected by BAT on 25 March 201 1 
with a mean count rate (integrated over 24 hours) of (0.0059 ± 0.0016) cts cm~ 2 s (too low to generate an on- 
board trigger). We note that the source was outside the BAT FOV for 2.86 hours before the first trigger. Between 
25 March and 31 March the BAT data show multiple flares peaking at up to 0.09 cts cm -2 s _1 (about 22 mCrabs). 
During the brightest parts of the outburst, the source had power law photon indices ranging from 1.3 to 1.8. The 
average 15-150 keV flux in the time interval 12:57:45 UT on 28 March 2011 to 05:30 UT on 30 March 2011 was 
about (1.0 ±0.2) x 10~ 9 erg cm~ 2 s _1 . The source flux then dropped dramatically, with an average count rate of 
0.002 ± 0.0005 cts cm" 2 s" 1 between 2 April 201 1 and 12 April 201 1 (see Supplementary FigureQ}- 

For most of the BAT data the spectral slope cannot be determined, making it difficult to determine an accurate 
flux based on the count rates. However if we assume a spectral slope, we can obtain estimates of the 15-150 keV 
flux. Supplementary Table |2]provides flux conversion factors for average BAT count rates at 3 epochs, for 2 observed 
spectral slopes. 

XRT data Observations of Sw J1644+57 with the Swift X-ray Telescope^ (XRT, 0.3-10 keV) began at 13:20:52 UT 
on 28 March 201 1. XRT observed Sw J1644+57 daily for between 1 1 ks and 28 ks per day. We report here on the first 
23 days after the trigger, during which the mean time spent on-target was 15.2 ks per day. As Swift is in a low-Earth 
orbit, observations are broken into small snapshots, typically 20-30 min long per ~ 96 min Swift orbit. 

For the first 14.5 days the XRT was in 'Auto State", in which it autonomously selects the appropriate observing 
mode based upon the brightness of the observed source. This meant that as the source varied in brightness, the XRT 
collected data in a combination of Windowed Timing (WT) and Photon Counting (PC) modes. PC mode is the standard 



Supplementary Table 1 | Swift BAT trigger details for Sw J1 644+57. 



Trigger Number Date Time Trigger Duration Intensity 



UT (s) (counts s ) 

450158 28 Mar 12:57:45 1208 6T 

450161 28 Mar 13:40:41 64 19.4 
Threshold set to so BAT would trigger on it again 

450257 29 Mar 18:26:25 320 15.6 

450258 29 Mar 19:57:45 64 38.2 
Triggers from this source disabled 



1 



Supplementary Table 2 | Flux conversion for observed BAT count rates for different spectral slopes 



Time Period 


BAT count rate 

(cts s~ ! cm -2 ) 


Photon Index 
(assumed) 


Estimated Flux 

(erg cm" 2 s~\ 15-150 keV) 


Pre-Trigger 


0.0059 


1.8 


5.1 x 10~ 1U 




0.0059 


1.3 


6.1 x 10~ 10 


Peak rate 


0.0900 


1.8 


7.8 x 10~ 9 




0.0900 


1.3 


9.3 x 10~ 9 


Late time 


0.0020 


1.8 


1.7 x 10~ 10 




0.0020 


1.3 


2.1 x 10- 10 



Supplementary Table 3 | Swift UVOT UV upper limits. 





Tq + — 2 days 




To + 2 — 15 days 




Filter 


Exposure (s) 


UL(3a) 


Exposure (s) 


UL(3d) 


white 


18359 


>23.6 


11221 


>23.5 


V 


2541 


>21.1 


no observations 




b 


1493 


>21.5 


no observations 




Ll 


31593 


>23.1 


37032 


>23.2 


wl 


3422 


>21.8 


82742 


>23.5 


m2 


3048 


>21.6 


38324 


>23.1 


w2 


2320 


>21.8 


24763 


>23.1 



XRT imaging mode, with a timing resolution of 2.5 s; if a source is brighter than ~ 0.5 XRT count s in this mode, 
pile-up correction must be used 3 1 . WT mode is a fast timing mode ( 1 . 8 ms time resolution) that collects 1 -dimensional 
image data. The fast readout of WT mode means that pile-up is avoided for sources with a count rate below ~ 100 
counts s _1 . 

The switch points for Auto State are tuned specifically for GRB observations, i.e. a fading X-ray light curve. For 
sources of moderate brightness (1-5 XRT counts s ), data will often be collected in PC mode. For this reason we 
changed our observing mode to WT for data collected from 12 April 201 1 onwards, to avoid issues relating to pile-up. 

X-ray spectra and light curves were produced utilizing the methods described by Evans et al. 32 All PC mode data 
were corrected for pile-up by removing events from the core of the PSF, and utilizing the wings of the PSF for both 
spectral fitting and light curve generation 31 . XRT spectra were corrected for the effects of charge traps that have 
developed due to radiation damage to the CCD over the Swift mission lifetime. 

The XRT count rate light curve for Sw J 1644+57 is shown in Supplementary Figure [2] Total exposure time for the 
first 23.4 days is 3.56 x 10 5 s, with a live time fraction of 17.6%. 

We produced light curves in two energy bands (soft and hard) in order to examine spectral variations with time and 
flux. We observe a strong anti-correlation in WT mode data between the spectral hardness and flux, with softer spectra 
as the source intensity decreases (see Supplementary Figure[3]l. There is a similar correlation in PC mode data, but the 
scatter is larger and the correlation is weaker. 

UVOT data The Swift UV-Optical Telescope^ 3 - (UVOT) observed Sw J 1644+57 in all optical and ultraviolet filters 
immediately following each of the four Swift BAT triggers, switching to mainly the ultraviolet filters for deep imaging 
after 1 April 2011. Observing periods were the same as for XRT. No persistent source has been detected in any 
filter. Photometry has been performed on individual and summed images using the UVOT photometric system—. 
No correction has been applied for Galactic extinction due to the reddening of Eb-v = 0.024 in the direction of the 
transient—. Supplementary Table [3] provides the upper limits for two epochs: early X-ray variability (first two days 
after the Swift BAT trigger), and later times. 
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Supplementary Figure 1 | BAT (top) and XRT (bottom) light curves of Sw J1 644+57. Times are 
measured from the BAT trigger on 28 March 2011. Although the higher count rate in the XRT allows 
better time resolution, it is clear that the count rates in both instruments track each other, with numerous 
flares in the first few days after the BAT trigger. Note that Sw J1 644+57 was detected by BAT —4-5 days 
before BAT triggered and XRT observations began. Data gaps are caused by periods when the source 
was not being observed. 



XRT Timing analysis 

After applying the barycentric correction to the XRT event lists we searched for coherent or quasi-periodic signals over 
different energy ranges and time intervals by using Fourier techniques. Due to the presence of strong non-Poissonian 
noise, mainly introduced by the rapid variability of the source, particular care must be taken in evaluating the statistical 
significance of any candidate signal. Following the prescriptions of Israel & Stella 36 , no significant (periodic or quasi- 
periodic) signal was found. Conservative 3o~ upper limits on the pulsed fraction, defined as semi-amplitude of the 
sinusoid divided by the mean source count rate, were computed. For the data span in the 31 March 201 1 to 19 April 
201 1 interval (i.e. without the initial intense spikes in order to mitigate the low frequency noise in the power spectrum) 
and the energy range 0.3 — lOkeV (524288 trials between 1 x 10~ 6 and 0.2 Hz), the 3(7 upper limits are in the l%-3% 
range for periods shorter than 500s, and in the 5%-90% range for periods between —500s and 1000s. 

X-ray spectral analysis 

The XRT data were processed with standard procedures (XRTPIPELINE vO.12.4), filtering and screening criteria by 
using FTOOLS in the Heasoft package (v.6.9). We used the latest (pre-release) version of the XRT software task 
xrtcalcpi, which calculates Pulse Invariant (PI) event energies taking into account position dependent corrections 
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Supplementary Figure 2 | Swift XRT count rate light curve of SwJ1 644+57, through 14 April 2011. 

Following nearly 3 days of intense flaring with peak count rates of over 80 cts s _1 , the source decayed 
for several days to a count rate of about 0.7 cts s -1 , then rose to —2 cts s~\ where it has remained (with 
excursions) for over 9 days. From —13 days after the BAT trigger, the light curve can be described by a 
mean level of —1 cts s -1 , with episodes of dipping that last 1 —2 days where the count rate drops to as 
low as 0.2 cts s _1 . 



for energy losses incurred by charge traps that have developed on the CCD. These software and calibration files 
are only available to the XRT team at the time of this analysis, but will be released to the public by the time this 
paper is published. They are the result of several years of calibration efforts by the XRT team to correct for the 
effects of radiation damage to the XRT detector after > 6 years in orbit. They provide significant improvements 
in the spectral resolution and flux calibration compared with earlier versions. We used the latest pre-release spec- 
tral redistribution matrices (swxwt0to2s6_20010101v013.arf and swxwt0to2s6_20070901v013.rmf for WT data and 
swxpc0tol2s6_20010101v012.arf and swxpc0tol2s6_20070901v012.rmf for PC data). 

Intensity-selected spectra We produced XRT spectra in seven intensity ranges for data collected from March 28 
until April 7 2011. Events were accumulated with count rates in the following intervals < 0.5, 0.5 — 1, 1 —2.5, 
2.5 - 10 counts s _1 for the PC data, and < 1 1, 1 1 - 35, > 35 counts s _1 for WT data. We corrected the PC data for 
pile-up by determining the size of an exclusion region at the core of the PSF necessary to get agreement between the 
wings of the observed PSF and the nominal PSF 31 , and excluding from the analysis all the events that fell within that 
region. 

The spectra were fit with several models. In all cases the absorption consisted of two components, one fixed to 
the Galactic value (Nn = 1.7 x 10 20 cm -2 ), and one free at the redshift z = 0.35. The following models provided 
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Supplementary Figure 3 | Ratio of XRT count rates in the 0.3-1.3 keV and 0.3-10 keV bands, 
plotted against the overall XRT count rate. The band ratio is strongly correlated with teh overall 
count rate, demonstrating that the spectrum gets progressively harder as the flux increases. There is a 
discontinuity around 4 ct/s; data at higher count rates were collected during flares in the first few days of 
monitoring, whilst data at lower count rates were collected after the XRT observing mode was changed 
from Auto to WT 14.5 days after the BAT trigger. The trend for harder spectra at higher count rates is 
seen on both sides of this discontinuity. This figure includes only WT-mode data taken within 23 days 
of the first BAT trigger. Each data point averages all WT data collected within a single Swift orbit. Only 
data points with uncertainties smaller than 25% are shown. 



good fits to the data: i) a simple absorbed power-law model (fbabs*zwabs(powerlaw)); ii) a log-parabola model 
(tbabs*zwabs(power2)); Hi) a broken power-law model (tbabs*zwabs(bknpower)). iv) an absorbed power-law model 
plus a diskblackbody (tbabs*zwabs(powerlaw+diskbb)). The "power2" log-parabola model is defined as: 



A{E)=E { -- a+ P* 10 ^ 



(1) 



The spectra are comparably well-fit well by either a log-parabola, a broken power-law or power-law plus multi- 
temperature thermal model. Again, we observe a trend for softer spectta as the source intensity decreases (see Supple- 
mentary Figure |4j. 



Time-selected spectra Spectra were also accumulated from all available XRT data with the aim of obtaining at 
least 3000 counts per specttum, unless relevant intensity variations were observed, in which cases more spectra were 
extracted, to represent the different intensity states. Hence some spectra were also exttacted with fewer counts than 
3000. We obtained a total of 32 PC spectra and 98 WT spectra, and we fit them with the first three models described 
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Supplementary Figure 4 | Swift XRT intensity-resolved spectral fitting results utilizing an ab- 
sorbed power-law model. Flux is for the 0.3-10 keV band. The power law slope is strongly correlated 
with the count rate, confirming the band ratio results shown in Supplementary Figure[3j 



above. The simple power-law model is an adequate description of the data, as is the log-parabola (F-test probabilities 
~ 10~ 3 ). We find that the column density is variable in time with no evidence of increased column when the dips occur. 
The results of the temporally resolved log-parabola model fits are shown in Supplementary Figure |5] in particular, the 
harder when brighter trend is also observed (the average photon index during the first three days and on the 15th day 
is 1 .78 ± 0.02 and 1.91 ± 0.08, respectively). 

Flux-calibrated XRT light curves The results of the spectral analysis with the log-parabola fits were used to cal- 
culate energy correction factors (ECFs) for conversion from the XRT count rate light curve shown in Supplementary 
Figure |2] to the flux-calibrated light curve presented in Figure 2 of the main Letter. The absorption-corrected fluxes 
(corresponding to the observed fluxes in Supplementary Figure [5]l were used to calculate ECFs corresponding to each 
time-resolved spectrum. These ECFs were then interpolated to the time of each bin in the count rate light curve of 
Supplementary Figure [2] The interpolated ECFs were multiplied by the counts in each light curve bin to obtain the 
absorption-corrected flux light curve. 

The flux light curve was used to estimate the total X-ray energy produced by Sw J1644+57. The total measured 
unabsorbed fluence is 1.24 x 10~ 4 erg cm~ 2 . Correcting for the live-time fraction, this gives a total unabsorbed fluence 
for Sw J1644+57 of 7.1 x 10 -4 erg cm~ 2 in the observed 0.3-10.0 keV band. Each flux point in the light curve was 
also converted to a luminosity at the source, and these were summed to obtain an estimate of the total X-ray energy 
(corrected for live-time fraction) of 3 x 10 53 erg in the 0.4-13.5 keV rest frame energy band. Mean, median, peak and 
minimum flux and luminosity values for Sw J1644+57 are given in Supplementary Table|4] 
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Supplementary Table 4 | X-ray Rest Frame Flux and Luminosity of Sw J1 644+57. 



Observer frame Rest frame 

Flux Level Unabsorbed Flux (0.3-10.0 keV) Luminosity (0.4-13.5 keV) 



Max flux 


1.03 x 


10 


^3 

erg 


cm" 


1 s" 1 


4.30 


X 


10 48 


erg 


;s-' 


Mean flux 


1.35 x 


10 


erg 


cm" 


- 2fj -l 


5.63 


X 


10 47 


erg 


;s-' 


Median flux 


7.07 x 


io- 


10 er f 


; cm 


"2 s- 1 


2.96 


X 


10 47 


erg 




Min flux 


3.35 x 


io- 


12 

er£ 


; cm 


"2 s- 1 


1.40 


X 


10 45 


erg 





XRT spectra for SEDs We also extracted spectra for use in our SED during three epochs representing three intensity 
states: i) strictly simultaneous with the event-mode BAT spectrum; ii) a very low state (4.5-5. 1 days after the trigger); 
Hi) an intermediate state (6.5-9.4 days after the trigger). The XRT data for these states are plotted in Supplementary 
Figures EMU 

Simultaneous BAT+XRT spectral fits The results of fitting the BAT+XRT spectral data simultaneously for one 
of the early flares can be seen in Supplementary Figure |6] We find that for this time period the BAT and XRT 
spectral fit are consistent, fitting well a broken power-law model. The BAT 15 — 150 keV spectrum is a continuation 
of the XRT spectrum, and there is no need for additional spectral breaks or high energy components to fit the BAT 
spectrum. Unfortunately it is not possible to fit perform simultaneous fits at all epochs, due to source faintness and 
the limitations of the Survey data collection mode of the BAT during regular observations. The spectral parameters 
appear to be strongly dependent on the X-ray flux (see Supplementary Figures |3]-|4}. Therefore in order to test whether 
the BAT spectra were typically consistent with the XRT fits, we fit a broken power-law model to the XRT spectrum at 
various flux levels, and then extrapolated this model to estimate the expected count rate seen in the BAT energy band. 
We found that the BAT count rates are highly consistent with the extrapolated XRT spectral fits, suggesting that the 
consistency seen in the BAT+XRT spectral fit show in Supplementary Figure |6]is true also for later time data. 

Observations by other X-ray/y-ray observatories 

MAXI Following the BAT on-board trigger, the MAXI 31 team reported that the Gas Slit Camera 38 also detected 
Sw J 1644+57—. The source rises to a peak about a day after the first BAT trigger. MAXI continued to detect the 
source for the initial bright flaring part of the outburst. The MAXI one-day-averaged light curve for the period before 
and after the BAT trigger is shown in Supplementary Figure [7] 

Fermi In the GeV energy range, no significant y-ray emission has been seen by the Fermi Large Area Telescope 
(LAT) from the direction of Sw J1644+57. During the period of main activity (28 March to 4 April 201 1) Fermi LAT 
observations constrain the average y-ray flux to < 2.7 x 10" 11 erg cm" 2 s~'(100 MeV-10 GeV, 95% confidence upper 
limit). During the following week, the equivalent upper limit is 2.9 x 10" 11 erg cm" 2 s -1 . Daily upper limits during 
the X-ray outburst vary, but are typically a few times 10~ 10 erg cm" 2 s , For the LAT analysis the event class used 
was "DIFFUSE" and the version of the instrument response functions used was Pass 6, v3 (P6_V3_DIFFUSE). 

UV/Optical/NIR observations 

Sw J1644+57 has been observed by a large number of ground-based optical and radio telescopes. In general, the 
heavily extincted optical counterpart is not clearly detected in optical or UV bands, where the host galaxy dominates 
the light, but is detected strongly in the near-infrared. Spectroscopic observations identified emission lines at a redshift 
of 0.35 1 4Q i 41 . Here we report details of optical and NIR observations obtained by our team. 

R-band Observations /?-band imaging data were taken using a CCD camera on the Mt. Lemmon Optical Astron- 
omy Observatory (LOAO) 1-m telescop e 42 ] 43 in Arizona, US, and SNUCAM on the Maidanak Observatory 1.5-m 
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Supplementary Figure 5 | Si/Wft XRT time resolved spectral fitting utilizing an absorbed log- 
parabola model. Parameters a and j3 for the log-parabola model are as defined in Supplementary 
Equation Q] Nh is the redshift corrected intrinsic absorption in units of 10 22 crrr 2 , XRT flux is the ob- 
served flux in units of erg crrr 2 s _l (0.3-10 keV). 

telescope in Uzbekistan.. The data were taken during the nights of 29 March 2011 to 8 April 2011 at LOAO, and 
12 April at Maidanak. A dithered sequence of 300 s - 600 s exposures were taken, resulting in 20min - 1 hour total 
integration for each night, which gave 3o~ detection limits of R = 22.4 — 23.9 mag over a 3" diameter aperture. The 
photometric calibration is based on field calibration performed using four Landolt standard star fields (PG1633+099, 
SA110, SA107 and PG1657+078) by the Lulin One-meter Telescope operated by National Central University. The 
summary of our observations and results are in Supplementary Table [5] 



BOAO/KASINICS and TNG observations and data analysis We imaged the field of S w J 1 644+57 with the Korea 
Astronomy and Space science Institute (KASI) Near Infrared Camera System (KASINICS 44 ) on the 1.8m telescope 
at the Bohyunsan Optical Astronomy Observatory (BOAO) in Korea, and with the NICS camera— on the Italian 
3.6m Telescopio Nazionale Galileo (TNG), located in La Palma, Canary Islands. Near-Infrared J, H and K/Ks-band 
observations were carried out on nights 30 and 31 March 2011 and 1, 3, 12, 13 and 21 April 2011. All nights were 
clear, with seeing in the range 1.1"— 1 .3". The complete observing log is reported in Supplementary Table[6] 

Image reduction was carried out using the jitter pipeline data reduction, part of the ESO-Eclipse package. As- 
trometry was performed using the 2MASS (http://www.ipach.caltech.edu/2mass) catalogue. Aperture photometry was 
made with the photom tool implemented in the GAIA package and the photometric calibration to the Ks band was 
done against the 2MASS catalogue. In order to minimize any systematic effect, we performed differential photometry 
with respect to a selection of local, isolated and unsaturated reference stars visible in the field of view. 
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Supplementary Figure 6 | Swift combined BAT/XRT "unfolded" spectra for Sw J1 644+57. The 

black points are XRT data and the red points are BAT data. The data are from the time interval 1 9:57:52.0 
- 20:04:32.0 UT on 29 March 2011, during one of the bright flares. This figure shows the spectrum 
incident on the instruments. The solid line is the fitted broken power-law model. The spectra are equally 
well-fit by a combination of a multi-temperature disk and a power law. 



The Swift XRT light-curve (Supplementary Figure [2) shows strong variability at essentially any epoch, and hints 
of an analogous behaviour in the NIR band are seen when considering the whole dataset. In particular, after the 
maximum NIR flux observed about 2 days after the first Swift trigger, global flux decreases are seen in the NIR data 
(see Supplementary Figure|8]l. About a week after the trigger the NIR flux seems to have entered a rather stable phase, 
which is still lasting more than three weeks after the burst. It seems likely that the late-time NIR data are due to the 
host galaxy contribution. 

By means of almost simultaneous z, J, H, K, and L observations (combining our data with those of Levan et 
al.— ) in the first days after the Swift trigger when the transient was brighter, we can use the SED of the transient to 
evaluate the presence of substantial rest-frame extinction. Modelling the SED with a power-law in frequency we find 
that solutions with a rather flat or even increasing spectrum are favored (v~ 1/3 ). The required extinction, assuming 
a Milky Way extinction curved, is in the range Eb-v ~ 1 — 3, with a strong covariance with the spectral index. The 
required extinction is indeed consistent with the observed constancy of the R-band flux. The R-band emission is due 
to the galaxy only and was already detected before the transient event 47 , while the transient contribution is depressed 
by local intrinsic absorption well below the host-galaxy brightness. 

Historical X-ray/y-ray Upper Limits 

We have examined archival databases from several X-ray and y-ray instruments that have observed the location of 
Sw J 1644+57, either through wide field surveying or in serendipitous pointings, to search for earlier detections of a 
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Supplementary Table 5 | Log of R-band observations of Sw J1 644+57. 



Start(UT) 



End(UT) 



Time(Mid Time) (UT) Filter Exp(s) Magnitude 



LOAO: 

2011-03-29T10:53:45 
2011-03-30T10:50:24 
201 1-04-05T1 1:15:26 
201 1-04-06T1 1:18:05 
2011-04-08T10:09:28 



2011-03-29T10:59:51 
201 1-03-30T1 1:47:47 
2011-04-05T12:02:33 
201 1-04-06T 12:22:00 
201 1-04-08T1 1:29:07 



2011-03-29T10:56:48 R 300x5 22.45 ± 0.49 

201 1-03-30T1 1:19:06 R 300x12 22.33 ± 0.22 

201 1-04-05T1 1:39:00 R 300x6 22.35 ± 0.25 

201 1-04-06T1 1:50:03 R 300x12 22.20 ±0.31 

2011-04-08T10:49:17 R 300x14 22.73 ± 0.26 



Maidanak Obervatory: 
2011-04-12T22:49:39 



201 1-04-12T23: 10:21 201 l-04-12T22:59:59 R 



600x2 22.49 ±0.12 



Supplementary Table 6 | Near Infrared data taken utilizing the TNG telescope and the BOAO telescope. 



Time of obs (UT) 


8t (days) 


Telescope 


Inst 


Filter 


Exp (s) 


Magnitude 


2011-03-30.650910 


2.11081 


BOAO 


KASINICS 


K s 


120x43 


16.78 ± 0.07 


2011-03-31.23164 


2.69558 


TNG 


NICS 


K 


20x3x40 


16.94 ± 0.04 


2011-03-31.636350 


3.09625 


BOAO 


KASINICS 


K s 


90x32 


17.07 ± 0.08 


2011-03-31.67763 


3.13753 


BOAO 


KASINICS 


H 


60x48 


18.44 ±0.09 


2011-04-02.13674 


4.59748 


TNG 


NICS 


J 


60x1x40 


20.10 ±0.12 


2011-04-02.17847 


4.63973 


TNG 


NICS 


H 


20x3x40 


19.12 ±0.09 


2011-04-02.22289 


4.68682 


TNG 


NICS 


K 


20x3x40 


17.77 ±0.10 


2011-04-03.648840 


6.10874 


BOAO 


KASINICS 


K s 


60x64 


18.09 ±0.11 


2011-04-03.689630 


6.14953 


BOAO 


KASINICS 


H 


60x96 


19.06 ±0.11 


2011-04-13.12765 


15.59158 


TNG 


NICS 


K 


20x3x20 


18.24 ±0.14 


2011-04-13.98145 


16.44759 


TNG 


NICS 


K 


20x3x40 


18.34 ±0.10 


2011-04-21.927812 


24.40270 


TNG 


NICS 


K 


20x3x30 


18.22 ±0.10 
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Supplementary Figure 7 | MAXI one-day-averaged light curve covering the —20 days before and 
after the detection of Sw J1 644+57 by BAT. 



source from this direction at a variety of wavelengths. 

BAT Historical upper limits The BAT instrument monitors about 80% of the sky daily. Examining the BAT Survey 
data at the location of Sw J1644+57, we find that the hard X-ray source never exceeded the 6a significance level when 
it was in the BAT field of view between February 2005 and March 2011 (count rate and flux limits vary with each 
observation), and was not found in any of the BAT survey studie s 48 ' 49 down to a limiting flux of 10~ n erg cm~ 2 s _1 
(15-150 keV), about two orders of magnitude below the average flux for the first 3 days after the BAT trigger. 

The BAT data are also monitored daily for transient events. We searched the entire BAT data set from 12 February 
2005 to 28 March 201 1, searching for previous detections of Sw J1644+57 on three timescales. The source is not de- 
tected on 16-day timescales, with a 3a upper limit of 0.001 1 cts cm~ 2 s _1 (about 5 mCrab). On one 1-day timescales, 
we obtain marginal detections exceeding 3a on 31 March 2009, 14 September 2009, and 14 March 201 1, all at about 
0.0036 ±0.001 1 cts cm~ 2 s _1 ; this rate is consistent with the expected rate of false 3a fluctuations in a search of 1866 
independent measurements, and we therefore do not consider these detections to be significant. With the exception 
of these three days, the 1-day 3 a upper limits vary depending on the total exposure of the source, but are between 
0.0010 and 0.0015 cts cm~ 2 s _1 for 90% of the observations. On shorter timescales, we can state that the source never 
exceeded the 6a significance level on any timescale between 64 s and 1200 s when it was in the BAT field of view 
(count rate and flux limits vary with each observation) until it triggered the BAT in a 1200 second image trigger on 28 
March 2011. 

ROSAT Historical Upper Limits Sw J1644+57 was serendipitously observed in a 6.5 ks ROSAT PSPC-B obser- 
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Supplementary Figure 8 | Optical/NIR light curves in the R (blue), J (green), H (brown) and Ks 
(red) bands of Sw J1644+57. No variability has been detected in the R band while in the NIR a decay 
trend roughly consistent with X-ray data is visible. Error bars are 1 a. 



vation made on 3 April 1992 covering the energy range 0.1 — 2.35 keV. The source is ~ 42' off axis and is affected 
by shadowing from one of the PSPC windows support ribs, giving an effective exposure time of 3.3 ks. The 3a 
upper limit 50 for the PSPC count rate is 8.3 x 10~ 3 cts s , corresponding to an observed 0.3 — 10 keV flux limit of 
< 2.9 x 10~ 13 erg cm~ 2 s _1 (observed) for any of the spectral models that fit the XRT data. We note that no source is 
detected at even the 1.5(7 level. This limit is more than an order of magnitude below the faintest portion of the XRT 
light curve, and four orders of magnitude below the brightest flares. Upper limits on the flux for different spectral 
models for Sw J 1644+57 are given in Supplementary Table [7] 

The ROSAT All-Sky Survey also covered this field with a relatively long exposure due to its high ecliptic latitude, 
collecting 940 s of data between 11 July 1990 and 13 August 1991. These data provide a 3(7 PSPC count rate upper 
limit™' 



'of 7.3 x 10~ 3 cts s" 1 , 



with a corresponding 0.3-10 keV flux limit of < 2.5 x 10 13 ergs cm 2 s 1 . 



XMM-Newton Historical Upper Limits Although there were no pointed observations of the Sw J1644+57 field 
by XMM-Newton prior to the outburst, between 26 August 2001 and 25 March 201 1 XMM-Newton slewed over this 
object on one occasion at 01:42:32 UT on 27 July 2005. The EPIC-pn routinely takes data during such slews 51 . 
No detection was made during the 3.9 s exposure made with the medium filter, resulting in a 2o~ upper limit of < 
0.79 cts s _1 (0.2-12 keV), calculated following the Bayesian approach of Kraft, Burrows & Nousek 50 . Unabsorbed 
0.2—12 keV flux limits were derived for the following three spectral models: a power law spectrum of photon index 
r = 2.8 absorbed by a Hydrogen equivalent column density of Nh = 7.7 x 10 21 cm~ 2 , a broken power law with 
T\ =4.98,T2 = 2.48,£ , f,, eflJ t = 1.96 keV, andAfy/ = 1.30 x 10 22 cm~ 2 , and a disk blackbody plus power law model with 
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Supplementary Table 7 | 3a upper limits based on ROSAT data for various spectral models 



All fluxes are 3(7 upper limits, erg cm 2 s '(0.3-10 keV) 

3 April 1992 Pointed Observation RASS Observation 

Model Observed Flux Unabsorbed Flux Observed Flux Unabsorbed Flux 



Simple absorbed power law model (Nh = 7.73 x 10 21 cm 2 ,T = 2.80): 

<2.88xl0~ 13 <1.21xl0~ 12 <2.52xl0~ 13 < 1.06x10" 



Broken power law model (A^ H = 1.30 x 10 22 cm~ 2 ,ri = 4.98,E break = 1.96 keV, F 2 = 2.48): 
<2.88xl0~ 13 <1.83xl0- n < 2.52 x 10~ 13 < 1.59 x 10~ n 



Disk black body + power law model (N H = 1-24 x 10 22 cm" 2 ,r in = 175eV,T = 2.40): 
<2.90xl0~ 13 <4.24xl0~ 12 



Ti„ = 175 eV, T = 2.40, and N H = 1 .24 x 10 22 cm" 2 ; these flux limits were < 2.9 x 10~ 12 erg cm" 2 s _1 , < 3.2 X 10~ 12 
erg cm -2 s -1 , and < 3.2 x 10~ 12 erg cm~ 2 s _1 , respectively. 

MAXI Historical Upper Limits To investigate whether MAXI detected the S w J 1 644+57 before the current outburst, 
we have analysed historical MAXI GSC data for the period of 17 August 2009 to 1 March 201 1 (see Supplementary 
Figure [9). No detection of Sw J1644+57 is found in those data, giving a 90% confidence level upper limit of ~ 
1.1 x 10~ n erg cm~ 2 s _1 (2-20 keV). In addition to this, utilizing data that has had additional cleaning applied, but 
covers a narrower time period ( between 1 September 2009 and 3 1 March 2010), we obtain a deeper 90% confidence 
level upper limit of 2.7 x 10 -12 erg cm~ 2 s (4-20 keV). 

Fermi Historical Upper Limits The analysis of the full 32-month Fermi/LAI dataset does not reveal any new 
source within 3 degrees of the position of the transient. This places 95% confidence upper limits of 1.7 x 10~ 8 
photons cnT 2 s -1 (100 MeV-10 GeV) and 1.5 x 10~ 10 photons cirT 2 s -1 (1 GeV - 300 GeV) on any persistent prior 
emission. Given the variable nature of this source, a search for emission was performed on shorter timescales of 2 
days and 5 days over the Fermi mission lifetime. No significant variation from the mean background level is observed 
in the light curve (see Supplementary Figure [Tot. 

Supplementary Discussion 

Constraints on the Mass of the Central Black Hole 

An upper limit to the mass of the central black hole in the galaxy where the transient was observed can be ob- 
tained through the black hole mass vs host bulge luminosity relation 52 . Using the measured B and H magnitudes of 
the host galaxy from published reports 53,54 and from our own analysis of Hubble Space Telescope WFC3 observa- 
tions, and taking a luminosity distance of of = 1.88 Gpc, we infer the following luminosities relative to the Sun: 
log(LB/i©,B) = 9.20 and log(Lft/L ( :),n) = 9.58 (the B-band luminosity includes a Galactic extinction correction of 
Aq = 0.08 mag). These luminosity estimates lead to a black hole mass of Mbh ~ 2 x 10 7 M Q , with a systematic un- 
certainty of a factor of 2 resulting from the scatter in the empirical relation between the host luminosity and black 
hole mass. Moreover, this empirical formula relates the luminosity of the spheroid of the host galaxy (i.e, its bulge) to 
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Supplementary Figure 9 | MAXI light curve of Sw J1 644+57 for the period August 17 2009 to 1 
March 2011. Note that the count rates seen in this light curve are consistent with a non-detection of 
Sw J1 644+57 in this period 

the black hole mass. However, the spectrum of the host galaxy shows the spectroscopic signature of star formation 55 , 
which suggests that it is a spiral. Since the magnitudes we have used refer to the entire galaxy, not just its bulge, the 
above value is an upper limit to the black hole mass. 

An independent constraint on the black hole mass can be obtained from the minimum X-ray variability time scale, 
which is observed to be 8t m m ~ 100 s (3(7 value). This time scale constrains the size of the central engine since, for 
Schwarzschild black hole, the minimum variability time scale is in its rest frame is 5f m i n ~ r$/c ~ 100 s, where 
M 6 =M bh /lO 6 M . Thus, 

M bh ~7.4xl0 6 (j|^) M . (2) 

Much larger or much smaller masses are unlikely for the following reasons. For larger masses, a variability time scale 
much shorter than the above constraint was never observed in the X-ray band in other systems containing supermassive 
black holes (e.g. AGNsfl For much smaller masses (e.g. stellar-mass black holes), one would expect to see much 
shorter variability time scales in the light curves. 

Some nearby low-luminosity GRBs with a stellar-mass central engine (e.g. GRB 060218) do have a smooth light 
curve without a noticeable fast variability component. However, GRB 060218 has only one peak. It would be contrived 
to have a smaller black hole smear up all the shorter variability time scales but only keep the 100 s time scale. We 
conclude that the black hole mass is likely M bh ~ (10 6 - 2 x 1O 7 )M (1 < M 6 < 20). 

'More rapid TeV variability has been observed in some blazars, which may be attributed to additional local Lorentz boost within a jet (e.g. jet 
in the jet)2£. 



14 



Supplementary Figure 10 | Fermi 95% confidence upper limit light curve of Sw J1 644+57 for the 
period August 17 2009 to 1 March 2011. 

Energetics and Their Implications 

After the initial flare, the source settled down to an X-ray flux of Fx ~ 10~ 9 erg cm~ 2 s -1 . The X-ray to y-ray 
luminosity is 

L r .x = 4.2 x 10 46 erg s" 1 ( F%X _ 2 ) , (3) 

\ 10 1U erg cm L s 1 / 

assuming the above luminosity distance. This luminosity is well above the Eddington limit for a 10 6 M black hole 
(^Edd = 1.3 x 10 44 M(, erg s _1 with Mbh = 10 6 Mg M ), suggesting that the emission originates from a relativistic jet 
rather than an accretion disk. 

The above luminosity, along with the coincidence of the event with the center of the host galaxy 54 , lead us to 
consider the onset of accretion-powered nuclear activity in the host galaxy. This could be fueled either by the tidal 
disruption of a star by the central black 57 hole or by the supply of gas from another source. 

To estimate the mass accretion rate and the total mass accreted, let us assume that the accretion rate is M = fM^dd, 
where MEdd = ^Edd/ (fjc 2 ), and rj is the usual accretion efficiency. We should allow the factor / to be larger than 
unity, since the accretion rate, especially for the early phases of the tidal disruption event, can be super-Eddington. 
Therefore 

M = /M Edd -/^ = 1.4xl0 23 M/gs- 1 (4) 
giving a total mass, accreted until now, of the order 

M = AtM ~2x 1CT 4 — — — — M Q (5) 
7] 1 month 

leading to a few per cent of a solar mass if f/t] ~ 10 2 . In fact, the measured radiative energy in the 0.4-13.5 keV 
rest frame band for the first 23.4 days of activity, 3 x 10 53 erg, corresponds to 0.016(?7 /O.1)M , assuming isotropic 
emission. 

If the source of fuel is the tidal disruption of a star by the black hole, then we expect an initially super-Eddington 
flare from the accretion of bound debris, followed by a steady decline that traces the rate of return of post-disruption 
debris, (M °< it — fo)~ 5,/3 ) 58 i 59 i 6Q . The monotonic decline phase could lag the initial flare by as much as 1-2 months 60 , 
though the X-ray light curve (Figure 1 in the Letter) suggests that the decline phase may have begun within a few days 
of the flare onset for Sw J1644+57. Therefore, long-term monitoring of the source in search of this specific rate of 
decline will be able to test the tidal disruption scenario. 
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Modelling the Spectral Energy Distribution of the Emerging Jet 
Available Data 

Radio Observations - The 8.4 GHz flux density was observed on 1 April 6 !, The 13.5 GHz, 15 GHz and 100 GHz 
(3 mm) fluxes were also reported in GCN Circulars 62 - 63 - 64 . Although important, these data will not be used to 
for the modelling of the X-ray and IR emitting region, since the radio emission must come from a larger region, 
not to be self-absorbed. Additional radio observations and their interpretation are treated by Zauderer et al.— . 

IR and Optical Observations - In the Supplementary Methods we present the near-infrared observations performed 
with the TNG and BOAO telescopes. Of these data, we have selected for the modelling those corresponding to 
the peak of the IR light curve (approximately 2 days after trigger), at the first minimum of the X-ray light curve 
(~4.5 days), and at an intermediate state about a week after the trigger. The large column density derived from 
the X-ray fitting, of the order of Nu ~ 10 22 cm~ 2 , suggests a large value of the extinction, even if the dust to 
gas ratio is not known. We de-redden our data with a rest frame Ay = 4.5, as discussed in the Supplementary 
Methods, approximately as in Bloom et al.— . 

X-Ray Observations - We have chosen to show the spectra for the peak (3 1 hours from trigger), including BAT data, 
for the minimum after 4.5 days, and for the intermediate persistent flux state that began a few days after the BAT 
trigger. 

y-ray Observations - We use the upper limit in the FermifLXT band reported in the Supplementary Methods. 

TeV Observations - The VERITAS Collaboration 67 set an upper limit on VHE emission during the early "flaring" 
period March 29-March 31 and the later "quiescent" period, April 1-April 15. The two upper limits are reported 
in the SED shown in Supplementary Figure [TTIjOl 

Assumptions and Observational Constraints on the Model 

We assume that the spectrum is produced by the synchrotron and inverse Compton (IC) mechanisms. In general, the 
seed photons for the inverse Compton process can be produced by the synchrotron process itself (Synchrotron Self- 
Compton, SSC) or can be produced outside the jet (External Compton, EC). A new-born blazar has had no time to 
build up a broad line region (i.e., to populate the surroundings of the jet with external photons), but some contribution 
to the EC process can come from photons produced in the accretion disk. The following is a summary of information 
we have used to guide the development of the model: 

• a limit on the size of the X-ray emitting region from the variability timescale, t v ~ 100 s; 

• the isotropic luminosity L bs; at the peak of the flare is L b s = 10 48 L syni 48 erg s _1 ; there is variability by a factor 
> 100; 

• some hints on the peak energy: Fermi and VERITAS did not detect the source; 

• the slope of the line connecting the IR to the X-rays is hard, possibly as hard as v 1 / 3 ; 

• it is the first time we have ever observed a flare of this kind, therefore the phenomenon is relatively rare. 

• the upper limit on the black hole mass derived above is 2 x 10 7 M and a mass compatible with the observed 
variability is — 7 x 10 6 M©. 

For blazar jets, one can find a unique solution for modelling their spectral energy distribution (SED) when the 
emission is SSC, and when we know the peak frequencies of the synchrotron and Compton spectra, the slopes before 
and after the peaks, the peak flux levels, and the variability timescale. Here we do not know anything about the 
Compton component, and we can only guess the peak frequency (and flux) of the synchrotron component. Therefore 
we need additional assumptions in order to find a reasonable fit. 
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One possibility is to assume that: (z) the source is magnetically dominated (note that this is not the case forblazars, 
whose jet powers are dominated by the kinetic energy of particles), and (if) the self Compton flux must be absorbed 
by the 7-7 — > e ± process. The first assumption comes from the very hard optical to X-ray spectrum, requiring a low 
energy cut-off in the particle distribution, and thus indicates a particle-starved source, while the second requirement 
comes from the Fermi and VERITAS upper limits. 

The size - We assume for simplicity a spherical source of radius R at some distance d from the black hole. We define 
l/[r(l — j3cosy/v)] as the Doppler factor, where y/ v is the viewing angle. Note that forblazars one assumes 
a mono-directional velocity (unlike in GRBs). The size is constrained by: 

R<ctv— (6) 
1 +z 

Magnetic field - We assume that the bulk of the luminosity is produced by synchrotron radiation. We further assume 
that the jet is magnetically dominated. This is because of the very hard slope between the optical and the X- 
rays, implying that there are no low energy electrons in the source. If there were, then they would severely 
overproduce the optical flux. One possibility is that the electrons are continuously reaccelerated/heated, in such 
a way that there are no electrons below, say, 7 ~ 10 2 - 10 3 . So we need a "thermal bath", but with relatively few 
electrons inside. Otherwise the mean energy of the particles will be small. As a consequence, the jet is particle- 
starved, and the jet power must be carried by the magnetic field. The jet is however capable of producing a lot 
of radiation: the isotropic luminosity we see is L c b s ~ 10 48 erg s -1 , but this is beamed. The real power spent by 
the jet to produce L D b s is 

P r ~^ (7) 

up to factor of order unity (see e.g. Ghisellini et al. 2010, MNRAS, 402, 497). If the jet is particle-starved, this 
power must be given by the Poynting flux, and therefore we require: 



Together with Supplementary Equation |6]this gives the limit: 

Pair production - Since we are dealing with a compact source, with a large produced synchrotron luminosity, it is 
conceivable that the self-Compton emission would be relevant. But the upper limits by Fermi tell us that the 
observed 7-ray flux is much less than the synchrotron flux. One way to account for that is pair opacity. We may 
therefore require that the optical depth for pair production is larger than unity: 



Or U' 

T ~?£^™ R> 1 (10) 
5 m e c z 

where t/ s ' sc is the comoving synchrotron self-Compton radiation energy density, and t/ s ' sc /m e c 2 is (very approx- 
imately!) the number density of photons at threshold, where Oyy ~ Ot/5. We then have: 

5 AnR 2 cm e c 2 207T L' syn Rm & c 3 20n L' syn @ 5 ct v m e c 3 

Distance from the black hole - If we assume a conical jet of semi-aperture angle = 0.1 1 we have 

, R 10ct v 3> 

e < el^Y+zj (12) 

Since we continue not to detect the source in 7-rays, while the X-rays vary widely, it is likely that we are 
observing some sort of a standing shock, or (magnetically dominated) internal shocks. The emission cannot be 
due to a single traveling blob. A single blob would travel rapidly, expand, and die in a short time. 
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Some numbers - From Supplementary Equation[TT]we can set a limit on 3i: 



9< nA f^K^c\ (13) 
y r v,2 ^syn / 

where f v = 10 2 f v ,2 s. Using this limit in Supplementary Equation|9] we derive a lower limit on B: 

r 1/10 2/i 
"syn,48'v,2 



B > 620 7W-^(# !1 ) Gauss (14) 



The size should be: 



R<2.5x lO'^L'^g cm (15) 



The distance from the black hole, in units of the Schwarzschild radius is 

A < 833 ^i« d6) 

Rs B-iM 6 
Modelling Procedure and Related Considerations 

We use the model described in detail in Ghisellini et al. 68 . The emitting region is assumed spherical, of size R^iob, 
moving with a bulk Lorentz factor T and located at a distance d from the black hole of mass M. The bolometric 
luminosity of the accretion disk is Ldisk- 

The particle energy distribution N(y) [cm -3 ] is calculated by solving the continuity equation where particle in- 
jection, radiative cooling and pair production (via the 7-7 — > e ± process) are taken into account. The created pairs 
contribute to the emission. The injection function Q(y) [cm~ 3 s _1 ] is assumed to be a smoothly joined broken power- 
law, with a slope Q(y) x Y~ Sl ar, d 7 -S2 below and above a break energy 7,: 

»' = e 'uW^ (17) 

In the specific application here, we assumed that electrons below a given 70 simply disappear. 

The total power injected into the source in the form of relativistic electrons is P[ = m e c 2 V J Q(y)ydy, where 
V = (4^/3)/?y G b is m e volume of the emitting region. 

The injection process lasts for a light crossing time ^biob/c, and we calculate N(y) at this time. This assumption 
comes from the fact that even if injection lasted longer, adiabatic losses caused by the expansion of the source (which is 
travelling while emitting) and the corresponding decrease of the magnetic field would make the observed flux decrease. 
Therefore the computed spectra correspond to the maximum of a flaring episode. 

Above and below the inner parts of the accretion disk there is an X-ray emitting corona of luminosity Lx (it is 
fixed at a level of 30% of La). Its spectrum is a power law of energy index ax ending with a exponential cutoff at 
E c =150keV. The specific energy density (i.e., as a function of frequency) of the disk and the corona are calculated 
in the comoving frame of the emitting blob, and used to properly calculate the resulting external inverse Compton 
spectrum. The internally produced synchrotron emission is used to calculate the synchrotron self Compton (SSC) 
flux. In this specific case we assume no Broad Line Region; the disk and the coronal radiation are negligible if the 
emitting region is at a large distance from the disk, but they become important producers of seed photons if the distance 
d of the emitting blob from the disk is of the same order of (or less than) the outer radius of the disk. 

Results and Discussion 

Supplementary Table [8] lists three sets of parameters adopted to generate three different models. In Supplementary 
Figures [TTlfl"3l we show the models resulting from these input parameters. Supplementary Table [9] lists the power 
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carried by the jet in the form of radiation (P r ), magnetic field (Pb), emitting electrons (P e , no cold electron component 
is assumed) and cold protons (P p , assuming one proton per emitting electron). All the powers are calculated as 

Pi = KR 2 bloh r 2 pcU! (18) 

where Uf is the energy density of the i component, as measured in the comoving frame. We note the following 
regarding Supplementary Table |9] 

• The power carried in the form of radiation, P T = ^biob^/^^rad' can be rewritten as [using U' aA = L' / [AnR\ Xoh c)\. 

r 2 r 2 1 

P r = L'—=L T ~L =■ (19) 

4 4^ 4 4^ 2 V ; 

where L is the total observed non-thermal luminosity (L' is in the comoving frame) and UL* is the radiation 
energy density produced by the jet (i.e. excluding the external components). The last equality assumes \jf v ~ 

i/r. 

• When calculating P e (the jet power in bulk motion of emitting electrons) we include their average energy, i.e. 
U{, = n e (y)m e c 2 . 

• For P p (the jet power in bulk motion of cold protons) we have assumed that there is one proton per emitting 
electron, i.e. electron-positron pairs are negligible. This is unimportant for the models shown in Supplementary 
FigureQT|and[T2] since in these cases the mean energy of the electrons is comparable to the rest mass of protons 
{(Y) m e ~ tn p ). However, for the model shown in Supplementary Figure [13] P p is the dominant form of power. 

• Pb is derived using the magnetic field found from the model fitting. 

In summary, we have studied the extreme cases of a jet whose power is largely dominated by the magnetic field 
(Models 1,2; Supplementary Figures ITTI and [TZt or by the kinetic energy of the matter (Model 3; Supplementary 
Figure Qj). 



A magnetic field-dominated jet 

Our rationale for this model is that since the accretion can be super-Eddington, the density of the accreting matter is 
large and can sustain a very large Z?-field, which launches the jet. In the dissipation region, the B-field is still large, 
possibly because it has not completed the acceleration of the matter. Dissipation can occur at the expense of the 
magnetic field, through reconnection. 

Models shown in Supplementary Figures QT| and [T2]rely on the SSC process (the magnetic energy density domi- 
nates over the radiation energy density of external seed photons), but adopt different black hole masses: 10 6 M and 
10 7 M Q , respectively. In the first model the disk luminosity is kept constant at the Eddington value, even if the jet 
luminosity and its power change, while in the second model we have assumed that the jet luminosity is tracking a 
fastly varying accretion luminosity. In this case, for the high state we have assumed a super-Eddington luminosity, and 
nearly Eddington for the low state. As can be seen, the disk luminosity can barely contribute to the soft X-ray flux. 

The dotted orange line in both models shows the flux emitted from a much larger region of the jet, producing the 
radio flux. 



A matter dominated jet 

For this model we assume that even if, at the start, the jet is dominated by the magnetic field, nevertheless the Poynting 
flux is able to accelerate the jet to its final F before the dissipation region. The jet is then dominated by the kinetic 
energy of matter, therefore dissipation is at the expense of the latter, as in all other blazars. The source of the seed 
photons is an accretion disk with an outer radius of > 300 Schwarzschild radii (it is 500 Rs for the model shown in 
Supplementary FigureQj}. This model is similar to models used for blazars, but with a particle energy distribution that 
does not extend up to such large energies (i.e. y max ~80-100), with most of the electrons at low energies. Therefore in 
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Supplementary Table 8 | Input parameters used to model the SED. Col. [1]: state; Col. [2]: size of the 
emitting region, in units of 10 15 cm; Col. [3]: black hole mass in solar masses; Col. [4]: Disk luminosity in 
units of 10 45 erg s _1 and (Ld/^Edd); Col. [5]: power injected in the blob calculated in the comoving frame, in 
units of 10 45 erg s -1 ; Col. [6]: magnetic field in Gauss; Col. [7]: bulk Lorentz factor; Col. [8]: viewing angle 
in degrees; Col. [9]: Doppler factor; Col. [10], [11] and [12]: minimum, break and maximum random Lorentz 
factors of the injected electrons; Col. [13] and [14]: slopes of the injected electron distribution [Q(y)\ below and 
above %\ Col. [15] f v = #biob(l + z)/(c@). The disk has an X-ray corona of luminosity L x = 0.3Ld- The spectral 
shape of the corona is assumed to be « v _1 exp(-/iv/150 keV) for the models in Supplementary Figure [TTIand 
Supplementary FigureQU while it is <* v~°- 7 exp(-/!v/150 keV) in Supplementary Fiaure[l2l 
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Model 2: Supplementary Figurell2l 
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Model 3: Supplementary Figure [13] 
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this case we have a sort of bulk-Compton process, in which relatively cold electrons scatter the seed photons coming 
from the outer radius of the disk. This model is similar to the one presented by Bloom et al.— , but in our case the 
IR and X-ray photons are produced in the same region. The large Lx/^ir ratio is due to the large ratio between the 
energy densities of the external radiation and the magnetic field. This requires a small magnetic field, and therefore a 
matter-dominated jet. 
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Supplementary Table 9 | Jet power in the form of radiation, Poynting flux, bulk motion of electrons and 
protons (assuming one proton per emitting electron). 



State 






logPr 


logPB 


logPe 


logPp 


Lo.3-io/Ltot 


Lis-iso/Ltot 


Model 1: 


Supplementary Fi| 


rureQj] 














High 






46.12 


45.65 


41.07 


41.05 


0.216 


0.318 


Low 






44.33 


45.65 


40.13 


40.55 


0.630 


1.20e-6 


Large 






40.82 


45.65 


41.02 


42.43 






Model 2: 


Supplementary Fi| 


!ure!12l 














High 






45.78 


46.65 


40.24 


40.34 


0.218 


0.354 


Low 






44.03 


44.78 


40.77 


40.95 


0.512 


8.5e-7 


Large 






41.30 


45.17 


41.31 


42.71 






Model 3: 


Supplementary Fij 


;ure[!3] 














High 






46.14 


42.47 


44.71 


47.74 


0.098 


0.175 


Low 






44.17 


41.64 


43.98 


47.19 


0.695 


0.038 
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Log v [Hz] 



Supplementary Figure 11 | Spectral Energy Distribution for Sw J1 644+57. The green data points show the 
transient brightness in the near-infrared (NIR) Ks band 48 hours after the trigger, while the cyan points show the 
J, H, and Ks band fluxes 4.5 days after trigger. The NIR flux has been dereddened with A v — 4.5. In the X-ray 
band we show the spectrum at the peak of the bright flares (31 hours after BAT trigger; XRT and BAT) and in the 
low flux state (~ 4.5 days after, XRT only), together with the spectrum in the intermediate persistent flux state 
that began a week after the BAT trigger (time integrated between day 6.5 and 9.5 from trigger, XRT only). The 
X-ray data have been corrected for absorption with a constant Ng = 2 x 10 22 crrr 2 . Upper limits from the Fermi 
LAT at 2 x 10 23 Hz and from VERITAS at 10 26 Hz 67 are also shown. The red curve shows the model discussed in 
the text, which is a blazar jet model 68 fit to our SED. The dominant emission mechanism is synchrotron radiation 
peaking in the X-ray band. On the low frequency side, the hard slope between the NIR and X-ray bands requires 
suppression of low-energy electrons, which would otherwise overproduce the NIR flux. On the high frequency 
side, the LAT and VERITAS upper limits require that the self-Compton component is suppressed by y-y pair 
production, without which the model would follow the dashed curve and would significantly overproduce the GeV 
and TeV emission. The model includes a disk/corona component from the accretion disk (black dotted curve), 
but the flux is dominated at all frequencies by the synchrotron component from the jet. The blue curve shows the 
corresponding model in the low X-ray flux state. The kink in the X-ray spectrum suggests a possible additional 
component may be required; it would have to be very narrow, and its origin is unclear. The radio fluxes come 
from a larger region of the jet (orange dotted line). See Supplementary Tableland Supplementary Table|9]for 
the model parameters of this fit. 22 




Log v [Hz] 



Supplementary Figure 12 | Data as in Supplementary Figure QT| Here the model is similar to that shown in 
Supplementary Figure Ql] but the black hole mass is larger (Mbh = 10 7 M ), and the jet luminosity is assumed 
to approximately track a changing disk luminosity, from 30xEddington in the high state to 0.6xEddington in the 
low state (dotted lines). The basic features of the jet are similar to the previous case, but the Poynting flux is not 
constant, instead it tracks M. 
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Log v [Hz] 

Supplementary Figure 13 | Data as in Supplementary FigureHU In this model the NIR flux is synchrotron 
radiation, and the X-rays are produced by the external Compton. The model assumes a matter dominated jet 
whose dissipating region is very close to the accretion disk (i.e. 120 Rs, the black hole mass is assumed to be 
1O 7 M ). The disk itself is assumed to extend out to —500 Rs- In this case the external part of the disk produces 
IR photons that can be efficiently scattered by the relativistic electrons in the jet. In order to avoid overproducing 
the NIR, the magnetic field is small, and the jet must be matter-dominated. The high state (red line) and the low 
state (blue line) differ in the amount of injected power in relativistic electrons and in the bulk Lorentz factor. The 
short-dashed gray line is the synchrotron flux in the high state, while the long-dashed line is the corresponding 
SSC component. 
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